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Description 

ITERATIVE CT RECONSTRUCTION 
METHOD USING MULTI-MODAL EDGE 

INFORMATION 

Background of the Invention 

[0001] The present disclosure relates generally to computed to- 
mography (CT) imaging and inspection systems and, more 
particularly, to a method for iterative CT reconstruction 
using multi-modal edge information. 

[0002] | n present date industrial inspection processes, different 
types of measurement systems are available such as com- 
puted tomography (CT), coordinate measuring machines 
(CMM), laser-based profilometry, etc. Each inspection 
modality has its own advantages and disadvantages asso- 
ciated therewith. Modalities such as CMM and laser-based 
profilometry can measure external surfaces with high ac- 
curacy, but cannot measure internal features unless the 
part is cut open. To date, CT is the most versatile of the 
measurement/inspection systems for revealing both the 



internal and external structures of industrial parts in a 
non-destructive manner. The potential industrial applica- 
tions of CT include reverse engineering, rapid prototyp- 
ing, casting simulation & validation, tire development, 
first article inspection, ceramic porosity inspection, pro- 
cess validation, parts qualification and defect detection, to 
name a few. However, the low inspection accuracy of CT 
inhibits widespread applications thereof. 
[0003] For example, in the area of reverse engineering, CT is a 
poor choice for capturing detailed external surface fea- 
tures, which can be crucial for capturing the design intent. 
The factors affecting CT accuracy in this regard include 
(among other aspects) beam-hardening, partial volume 
effect, scattering and off-focal radiation. Thus, in order to 
improve CT inspection accuracy, more effective methods 
are needed for removing the effects of these artifacts. In 
the area of CT image reconstruction, filtered backprojec- 
tion (FBP) is a common technique because of its fast com- 
putation and ease of implementation. However, because 
FBP oversimplifies the CT data acquisition into an ideal 
Radon transform (i.e., Fan Beam transform, cone beam 
transform or any other transform depending on the par- 
ticular acquisition geometry), the reconstructed image 



suffers from artifacts such as beam hardening and partial 
volume as discussed above, contributing to degraded im- 
age quality and inspection accuracy. Precorrection of the 
data prior to performing an FBP reconstruction is generally 
not capable of fully correcting for the artifacts due to 
these effects. 

[0004] Another type of CT image reconstruction methodology is 
what is referred to as iterative reconstruction. These tech- 
niques are based on different mathematical principles, 
such as the statistical approach of maximum likelihood, 
and the least squares approach, for example. Iterative 
methods allow for the incorporation of a dedicated for- 
ward model of the data acquisition. On the other hand, 
because the computation associated with iterative recon- 
struction is highly intensive, iterative methods are not yet 
widely used in CT. As such, in present inspection pro- 
cesses, CT images are typically reconstructed using the 
FBP method where both external and internal features are 
extracted from the CT images. In the event a more accu- 
rate measurement of internal features is desired, the part 
is generally cut open and inspected by CMM. Thereafter, 
the CT measurement is calibrated using the CMM mea- 
surement, so as to compensate for any possible bias. Un- 



fortunately, this procedure is both time-consuming and 
expensive. 

[0005] Accordingly, it is desirable to be able to provide an im- 
proved inspection method that accurately captures both 
the internal and external features of an object to be in- 
spected, but in a time efficient and inexpensive manner. 
Brief Description of the Invention 

[0006] The above discussed and other drawbacks and deficien- 
cies of the prior art are overcome or alleviated by an accu- 
rate computed tomography (CT) reconstruction method. In 
an exemplary embodiment, the method includes imple- 
menting an iterative image reconstruction process for CT 
metrology of an object, wherein the iterative reconstruc- 
tion process utilizes an accurate forward projection 
model. 

[0007] | n another aspect, a method for implementing an iterative 
reconstruction of a computed tomography (CT) image in- 
cludes, during each of a plurality of iterations, generating 
a reconstructed image. The reconstructed image is con- 
strained by utilizing prior outer edge information obtained 
from a modality in addition to CT, and then transformed 
to a projection domain so as to generate a calculated 
sinogram. A correction image is determined based on the 



calculated sinogram and a measured sinogram. 
[0008] | n s ti|| another aspect, a method for implementing an iter- 
ative reconstruction of a computed tomography (CT) im- 
age includes, during each of a plurality of iterations, gen- 
erating a reconstructed image. The reconstructed image is 
constrained utilizing prior outer edge information ob- 
tained from a modality in addition to CT and transformed 
to a projection domain utilizing so as to generate a calcu- 
lated sinogram. A correction image is determined based 
on the calculated sinogram and a measured sinogram, 
wherein the iterative reconstruction utilizes an accurate 
forward projection in determining the calculated sino- 
gram. 

Brief Description of the Drawings 

[0009] Referring to the exemplary drawings wherein like ele- 
ments are numbered alike in the several Figures: 

[0010] Figure 1 is a schematic block diagram illustrating a 
method for iterative CT reconstruction using a priori 
multi-modal edge information, in accordance with an em- 
bodiment of the invention; 

[0011] Figure 2 is a schematic block diagram illustrating a gen- 
eral iterative reconstruction process; 

[0012] Figure 3 is a schematic block diagram illustrating an itera- 



tive reconstruction scheme with accurate forward projec- 
tion, in accordance with a further aspect of the present in- 
vention; 

[0013] Figure 4 is a schematic block diagram illustrating an itera- 
tive reconstruction scheme with multi-modal constraint, 
in accordance with a further aspect of the present inven- 
tion; 

[0014] Figure 5 is a schematic block diagram illustrating an itera- 
tive reconstruction scheme with accurate forward projec- 
tion and multi-modal constraint as shown in Figure 1, in 
accordance with a further aspect of the present invention; 

[0015] Figure 6 is a schematic depiction illustrating the use of 

outer boundary information in the iterative reconstruction 
process of Figures 4 and 5; 

[0016] Figure 7 is an image produced using standard FBP recon- 
struction of simulated data, wherein beam hardening, 
partial volume and other artifacts are present; and 

[0017] Figure 8 is an image produced using the iterative FBP re- 
construction with accurate forward projection and multi- 
modal constraint. 
Detailed Description of the Invention 

[0018] Disclosed herein is a computed tomography (CT) recon- 
struction method having improved inspection accuracy, in 



which accurate forward projection is used in an iterative 
image reconstruction process applied to CT metrology of 
an object. Thereby, a more accurately calculated image 
projection (sinogram) is used in the iterative correction 
process. 

[0019] | n another aspect, external edge information from modal- 
ities such as coordinate measurement machine (CMM), 
micrometer and laser-based profilometry, (which are 
much more accurate than CT) or any other imaging device 
that gives some edge information about the inspected ob- 
ject, are incorporated into iterative CT image reconstruc- 
tion algorithms. The problems addressed by the present 
disclosure include the improvement of the overall accu- 
racy of the contours, the correction of artifacts and/or the 
improvement of the CT image quality, by using prior edge 
information. In addition to accurately modeling the 
physics of the acquisition, the invention embodiments 
herein may also utilize prior information (e.g., such as the 
edge information from a much more accurate measure- 
ment system) as constraints to correct artifacts, improve 
CT image quality, and/or improve the overall accuracy of 
the contours. The use of the prior edge information with 
iterative reconstruction is broadly applicable to the field of 



CT imaging. Furthermore, such prior information obtained 
from other modalities may be used in conjunction with 
accurate forward projection in the field of CT metrology. 

[0020] The particular reconstruction algorithms used may be ei- 
ther iterative filtered back projection (IFBP), or any other 
iterative algorithm such as those based on maximum a 
posteriori probability (MAP), maximum likelihood (ML), al- 
gebraic reconstruction technique (ART), entropy-based 
optimization, least squares (LS) or penalized weighted 
least squares (PWLS). In the embodiments utilizing prior 
edge information, the outer edge is constrained to satisfy 
the information obtained from the other modalities. In 
practice, this may be carried out during each iteration by 
adjusting the outer edge according to the information ob- 
tained by CMM or other laser-based profilometry, each 
time an intermediate reconstruction is obtained. 

[0021] Referring initially to Figure 1, there is shown a block dia- 
gram illustrating an overall process 100 for iterative CT 
reconstruction of an imaged object using a priori multi- 
modal edge information, in accordance with an embodi- 
ment of the invention. In particular, prior edge informa- 
tion 102 is obtained from an accurate measurement sys- 
tem such as a coordinate measuring machine (CMM) or 



laser-based profilometry, for example. This prior infor- 
mation 102 is used during the CT image reconstruction 
process, along with a measured sinogram 104 to remove 
artifacts or improve the edge position accuracy, and dur- 
ing the edge detection process to reduce the edge detec- 
tion error. As shown in Figure 1, an iterative image recon- 
struction scheme 106 receives the sinogram 104 and the 
prior edge information 102 as inputs thereto. Using both 
sinogram information and prior edge information, the ap- 
plication of the iterative image reconstruction results in an 
improved CT image, as shown in block 108. 

[0022] Moreover, the process 100 may also include post- 
processing steps that utilize the prior edge information 
102 and the improved CT image 108 for optimal thresh- 
old/adaptive edge detection, as shown at block 110. As a 
result, a two or three-dimensional point cloud 112 may 
be constructed for providing a description of the edges of 
the imaged object. 

[0023] | n order to appreciate the benefits provided by the present 
invention embodiments, reference is made to Figure 2 in 
which there is shown a schematic block diagram outlining 
a general iterative reconstruction process 200. In such a 
process, the image estimate is updated and transformed 



back and forth between the projection domain and the 
image domain. The particular manner in which transfor- 
mation between the domains is performed distinguishes 
different reconstruction methods. For instance, iterative 
FBP methods use FBP to reconstruct an image and they 
use a Radon or Fan Beam (or other, depending on the ge- 
ometry) transform to calculate a sinogram. 
[0024] M 0re specifically, a calculated sinogram 202 (i.e., a set of 
predicted projection data) is compared with the actual 
measured sinogram data 104, as shown at block 204. Ini- 
tially, the calculated sinogram 202 is a transformation 
from a predicted reconstructed image 206. In the first it- 
eration, the reconstructed image may be any appropriate 
default setting. As a result of the comparison of the cal- 
culated sinogram 202 and the measured sinogram 104, a 
sinogram error 208 is produced. The sinogram error 208 
is then transformed back to the image domain to produce 
a correction image 210. Thus, in a subsequent iteration, 
the correction image 210 is used along with the previous 
reconstructed image to create an updated reconstructed 
image 206 for the next iteration, as shown at block 212. 
Thereafter, the process is repeated for a desired number 
of iterations or until some other stopping criterion is met. 



[0025] As indicated previously, and with particular regard to cap- 
turing detailed external features, the reconstructed image 
206 suffers from artifacts due to physical effects such as 
beam hardening and partial volume. Therefore, as shown 
in Figure 3, an iterative reconstruction scheme 300 is dis- 
closed that includes an accurate forward projection model 
302 for use in CT metrology of an object. With the accu- 
rate forward projection (instead of using an ideal Radon 
transform), physical effects of the acquisition, such as a 
polychromatic x-ray spectrum, a finite detector point- 
spread function, a finite focal spot size, azimuthal blur, 
scattered radiation, measurement noise, or off-focal radi- 
ation, are each taken into account. In particular, an itera- 
tive maximum-likelihood polychromatic algorithm for CT 
(IMPACT) is used for preventing beam hardening effects. 
In the IMPACT algorithm, the continuous spectrum of the 
x-ray tube is modeled as a number of discrete energies. 
The energy dependence of the attenuation is taken into 
account by decomposing the linear attenuation coefficient 
into a photoelectric component and a Compton scatter 
component. The relative weight of these components is 
constrained based on prior material assumptions. Addi- 
tional information regarding the IMPACT algorithm as part 



of a forward projection model may be found in Bruno De- 
Man, Johan Nuyts, Patrick Dupont, Guy Marchal, and Paul 
Suetens, "An Iterative Maximum-Likelihood Polychromatic 
Algorithm for CT," IEEE Trans. Med. Imag., vol. 20, No. 10, 
pp. 999-1008, October, 2001, the contents of which are 
incorporated herein by reference in their entirety. 

[0026] | n addition, the iterative reconstruction algorithm may be 
adapted to take into account scatter. If a scatter estimate 
is available (from a measurement or a simulation), instead 
of subtracting it from the measurements, it can be used 
as an input of the reconstruction algorithm. Additional in- 
formation regarding scatter may be found in Bruno De- 
Man, "Iterative Reconstruction for Reduction of Metal Arti- 
facts in Computed Tomography," PhD Thesis, Katholieke 
Universiteit Leuven, May 21, 2001, the contents of which 
are incorporated herein by reference in their entirety. 

[0027] Furthermore, since finite beamwidth is one of the limiting 
factors in determining achievable spatial resolution of the 
scanner, it is also taken into account in the iterative re- 
construction algorithm by incorporation in the accurate 
forward projection, which would otherwise operate under 
the assumption that the beamwidth is infinitesimal. For 
instance, the accurate forward projection can model the 



detector point spread function, the finite detector aper- 
ture, the finite focal spot size, azimuthal blur, detector 
cross-talk, and off-focal radiation. Additional information 
regarding consideration of finite beamwidth in maximum 
likelihood algorithms may be found in Jolyon A. Browne, 
John M. Boone, and Timothy J. Holmes, "Maximum-Likeli- 
hood X-ray Computed-Tomography Finite-Beamwidth 
Consideration," Applied Optics, vol. 34, No. 23, pp. 
5199-5207, August 10, 1995, the contents of which are 
incorporated herein by reference in their entirety. 

[0028] Referring now to Figure 4, there is shown an alternative 
iterative reconstruction scheme 400 utilizing multi-modal 
constraint information 402, in accordance with a further 
aspect of the present invention. As is shown, the multi- 
modal constraint information 402 (e.g., external edge in- 
formation, such as obtained through CMM or laser-based 
measurement systems) is used to modify the current esti- 
mated reconstruction image, prior to calculating an up- 
dated sinogram 202 for comparison with the measured 
sinogram 104. In this exemplary embodiment, the multi- 
modal constraint information 402 is used to adjust the 
outer edge information in each iteration. 

[0029] Figure 5 illustrates still another iterative reconstruction 



scheme 500 in which both the multi-modal constraint in- 
formation 402 and the accurate forward projection 302 
model parameters are utilized in the sinogram calculation. 

[0030] Figure 6 is a schematic depiction illustrating the use of 

outer boundary information in the iterative reconstruction 
process of the embodiments of Figures 4 and 5. As is 
shown, the outer edges from an intermediate reconstruc- 
tion (iteration) 602 are combined with the outer edge data 
606 obtained from another modality to result in an up- 
dated image 608. In each successive iteration, outer edges 
from intermediate reconstruction are again adjusted by 
using the accurate outer edge information. 

[0031] Finally, Figures 7 and 8 illustrate a comparison (based on 
simulated data) between an image produced using stan- 
dard FBP reconstruction, and an image produced using 
accurate forward projection and multi-modal constraint, 
as described above. As can be seen from Figure 7, the 
conventional FBP process (with beam hardening and par- 
tial volume and other artifacts) results in an external edge 
RMS error of 1.04 mil and an internal edge RMS error of 
0.64 mil. On the other hand, Figure 8 is an image pro- 
duced using the iterative FBP reconstruction with accurate 
forward projection and multi-modal constraint. Compared 



with conventional FBP method in Figure 7, the present 
method further reduces the RMS error from 1.04 mil (FBP) 
to 0.06 mil for external features and 0.64 mil (FBP) to 
0.25 mil for internal features. 

[0032] As will be appreciated, the use of an accurate forward 

projection in each iteration removes image artifacts such 
as beam hardening and partial volume. Accurate forward 
projection also improves the image quality, and improves 
CT inspection accuracy. In addition, the use of external 
edge information in each iteration may also help the arti- 
fact correction procedure and improve the overall accu- 
racy of the contours. Another benefit is realized from the 
substantially increased convergence speed of the iterative 
reconstruction, by the use of outer edge constraints. Fur- 
thermore the embodiments described herein should not 
be construed as being limited to a specific type of CT ge- 
ometry. For example, the CT geometry used to perform 
the image reconstruction may be 2D or 3D, whereas the 
beam geometry may be parallel-beam, fan-beam, cone- 
beam, or any other exotic geometry types known to one 
skilled in the art. 

[0033] As will be also appreciated, the above described method 

embodiments may take the form of computer or controller 



implemented processes and apparatuses for practicing 
those processes. The disclosure can also be embodied in 
the form of computer program code containing instruc- 
tions embodied in tangible media, such as floppy 
diskettes, CD-ROMs, hard drives, or any other computer- 
readable storage medium, wherein, when the computer 
program code is loaded into and executed by a computer 
or controller, the computer becomes an apparatus for 
practicing the invention. The disclosure may also be em- 
bodied in the form of computer program code or signal, 
for example, whether stored in a storage medium, loaded 
into and/or executed by a computer or controller, or 
transmitted over some transmission medium, such as over 
electrical wiring or cabling, through fiber optics, or via 
electromagnetic radiation, wherein, when the computer 
program code is loaded into and executed by a computer, 
the computer becomes an apparatus for practicing the in- 
vention. When implemented on a general-purpose micro- 
processor, the computer program code segments config- 
ure the microprocessor to create specific logic circuits. 
[0034] while the invention has been described with reference to a 
preferred embodiment, it will be understood by those 
skilled in the art that various changes may be made and 



equivalents may be substituted for elements thereof with- 
out departing from the scope of the invention. In addition, 
many modifications may be made to adapt a particular 
situation or material to the teachings of the invention 
without departing from the essential scope thereof. 
Therefore, it is intended that the invention not be limited 
to the particular embodiment disclosed as the best mode 
contemplated for carrying out this invention, but that the 
invention will include all embodiments falling within the 
scope of the appended claims. 



